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ABSTRACT 

We describe a new method to measure the escape fraction / esc of ionizing radiation from distant 
star- forming galaxies using the afterglow spectra of long-duration 7-ray bursts (GRBs). Optical 
spectra of GRB afterglows allow us to evaluate the optical depth of the host ISM, according to the 
neutral hydrogen column density A (HI) observed along the sightlines toward the star- forming regions 
where the GRBs are found. Different from previous effort in searching for faint, transmitted Lyman 
continuum photons, our method is not subject to background subtraction uncertainties and does not 
require prior knowledge of either the spectral shape of the host galaxy population or the IGM Lyo 
forest absorption along these GRB sightlines. Because most GRBs occur in sub-L* galaxies, our study 
also offers the first constraint on / osc for distant low-mass galaxies that dominate the cosmic luminosity 
density. We have compiled a sample of 28 GRBs at redshift z > 2 for which the underlying AT(HI) in 
the host ISM are known. These GRBs together offer a statistical sampling of the integrated optical 
depth to ionizing photons along random sightlines from star-forming regions in the host galaxies, and 
allow us to estimate the mean escape fraction (/ CS c) averaged over different viewing angles. We find 
(/esc) = 0.02 ±0.02 and place a 95% c.l. upper limit (/ eS c) < 0.075 for these hosts. We discuss possible 
biases of our approach and implications of the result. Finally, we propose to extend this technique for 
measuring (/ CS c) at z ~ 0.2 using spectra of core-collapse supernovae. 

Subject headings: cosmology: observations — gamma-rays:bursts — galaxies:high-redshift — galaxies:ISM 



1. INTRODUCTION 

Observations of distant QSOs indicate that the inter- 
galactic medium (IGM) became fully ionized by redshift 
z ~ 6 (e.g. Fan, Carilli, & Keating 2006). While at z < 3 
QSOs are the dominant sources of the ultraviolet back- 
ground radiation (e.g. Haardt & Madau 1996), at higher 
redshifts where the number density of QSOs declines 
steeply toward earlier epochs (e.g. Willott et al. 2005; 
2005; Richards et al. 2006) additional ionizing sources 
are necessary. The spectral shape of the ultraviolet back- 
ground radiation inferred from intervening metal absorp- 
tion line studies (e.g. Haehnelt et al. 2001) and inter- 
galactic Hell absorption spectra (e.g. Shull et al. 2004; 
Reimers et al. 2005) suggest that young stars may pro- 
vide the dominant ionizing sources during early epochs. 

The escape fraction of ionizing radiation, / csc , speci- 
fies the fraction of stellar-origin ionizing photons (hv > 1 
Ryd) that escape the interstellar medium (ISM) of star- 
forming galaxies. Accurate measurements of / esc are im- 
portant for quantifying the relative contribution of ioniz- 
ing photons to the ultraviolet background radiation be- 
tween galaxies and AGN. For the Milky Way, estimates 
based on diffuse Ha emission of High Velocity Clouds 
(HVC) yield an upper limit of / osc < 6% (e.g. Bland- 
Hawthorn & Maloney 1999 and see Weiner et al. 2001 
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for a review). In the nearby universe, early observations 
of starburst galaxies place constraints at / csc < 6% (e.g. 
Heckman et al. 2001). Bergvall et al. (2006) have re- 
ported the first positive detection of Lyman continuum 
photons in the spectra of the starburst galaxy Haro 11 
using the Far-UV Space Explorer. Their analysis sug- 
gests an escape fraction in this galaxy of / csc = 1 — 10%. 
This finding has, however, been challenged by Grimes 
et al. (2007), who cannot confirm the detection of Ly- 
man continuum in the same data set. 

At higher redshift, a range of values are reported from 
/osc < 6% at z = 1.1 - 1.4 (Malkan et al. 2003; Siana 
et al. 2007), to between mean values of no more than 
8% (Giallongo et al. 2002; Fernandez-Soto et al. 2003; 
Shapley et al. 2006) and / osc ~ 13 — 38% (Inoue et al. 
2005) at z ~ 3, to / esc » 22% at (z) = 3.4 (Steidel et al. 
2001) 6 . In addition, direct detections of Lyman contin- 
uum photons are reported for two z ~ 3 galaxies by 
Shapley et al. (2006), implying / csc « 13 — 20% for these 
two sources. The large scatter in the reported / esc may 
imply a large variation in the optical depth across differ- 
ent lines of sight toward the inner regions of the galaxies 
(e.g. Gnedin, Kravtsov, & Chen 2007), but it also under- 
scores the challenges in detecting low-luminosity features 
in a background noise limited regime. 

In this Letter, we describe a novel approach for con- 
straining / osc from high-redshift star-forming galaxies. 
We estimate / esc based on the distribution of neutral hy- 
drogen column density iV(H I) observed in the afterglow 
spectra of long-duration 7-ray bursts (GRBs). Long- 

6 Note that many of the previous publications reported measure- 
ments for / esc f e i , which is defined as the ratio of escaped ionizing 
photons at 912 A to the observed flux density at 1500 A. Here we 
have converted these relative measurements to / esc at 1 Ryd based 
on their respective (E(B — V)) reported by the authors. 
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duration GRBs are believed to originate in the death of 
massive stars with M > 20 Mq (see Woosley & Bloom 
2006 for a recent review) , and are signposts of active star- 
forming regions in the ISM of their host galaxies (e.g. 
Bloom et al. 2002; Fruchter et al. 2006). Spectroscopic 
observations of the bright optical afterglows following the 
initial bursts have allowed us to measure the gas and dust 
content along the sightlines toward the GRBs, based on 
absorption features imprinted in the afterglow spectra 
(e.g. Vreeswijk et al. 2004; Chen et al. 2005; Prochaska 
et al. 2007a). These GRB sightlines together offer a sta- 
tistical sampling of the integrated optical depth to ioniz- 
ing photons along random directions in the host galaxies. 

The escape fraction / esc determined from the iV(HI) 
distribution along random sightlines in the ISM of GRB 
host galaxies does not require direct detection of Lyman 
continuum photons. It is not subject to systematic uncer- 
tainties due to background subtraction. The / CS c value is 
derived based on the total gas column observed in front of 
the star-forming region that hosts the GRB. It does not 
depend on the spectral shape of the ultraviolet radiation 
from the host galaxy or the stochastic uncertainties in the 
IGM Lya forest absorption along the lines of sight. Fi- 
nally, while some GRB host galaxies are reported to have 
high star formation rate, SFR > 100 M , (c.f. Berger 
et al. 2003; Le Floc'h et al. 2006), growing evidence in- 
dicates that the majority are sub-L* galaxies (e.g. Le 
Floc'h et al. 2003; Sollerman et al. 2005; Fruchter et al. 
2006). Our study therefore offers the first constraint on 
/osc for low-mass galaxies at z > 2. 

2. THE SAMPLE OF GRBS AT Z > 2 

To obtain an accurate estimate of the mean escape 
fraction of ionizing photons along GRB sightlines, we 
first compile a sample of GRBs that are confirmed at 
z > 2. We focus our analysis on bursts at zqrb ^ 2 
for two main reasons. First, it minimizes the ambiguity 
between an optically thin sightlinc and a low-redshift in- 
terloper. The z > 2 IGM will imprint its signature on an 
optical spectrum via the Lya forest and metal-line ab- 
sorption features. Therefore, an afterlow that exhibits a 
featurless spectrum is most likely at zqrb < 2. Second, 
measurements of iV(HI) from the absorption profiles of 
Lya and Lyman series allow us to directly evaluate the 
optical depth at the Lyman limit frequency ill- 

Table 1 lists 40 spectroscopically confirmed GRBs at 
zgrb ^ 2. In addition to the redshift of each source, we 
also list the isotropic equivalent energy release in 7-ray 
photons (£iso), the observed iV(HI) if available, a flag 
f Q to indicate whether the afterglow spectrum covers the 
rcdshifted Lya transition ('0' means no coverage and '1' 
means Lya coverage), and a flag fj to indicate whether 
metal-line features due to low ions such Si II or C II are 
present. Four of the GRBs do not have spectral cover- 
age of the Lya transition from the host, but the presence 
of low ions (Ji) indicate that the gas is consistent with 
being optically thick. To exhibit strong low-ion absorp- 
tion, an optically thin gas would need to have very high 
density (to maintain a non-negligible neutral fraction) 
and super-solar metallicity, both of which are very un- 
likely. Nine sources do not have published ./V(HI). In the 
subsequent analysis, we consider only those 28 sightlines 
with published iV(HI) values as our main sample, and 
assume that the remaining 13 sources with no available 




Fig. 1. — Cumulative distribution of neutral hydrogen column 
density T[< iV(HI)] observed in the host galaxies of long-duration 
GRBs at z > 2 (solid histogram). The shaded area shows the 
1-a uncertainties evaluated using a bootstrap re-sampling method 
that accounts for both 7V(HI) measurement uncertainties and sam- 
pling errors. The dash-dotted line represents the best-fit power-law 
model described in § 3. The dotted curve respresents the predicted 
distribution from Gnedin, Kravtsov, & Chen (2007). 

N(H I) measurements share the same distribution as the 
main sample. This is justified based on the similar zqrb 
and Ei so distributions between GRB sightlines with and 
without known iV(HI) measurements. 

Figure 1 presents the cumulative 7V(HI) distribution, 
J-[< N(H I)] from the main sample of 28 GRB host galax- 
ies, together with the l-er uncertainties determined based 
on a bootstrap re-sampling method. Specifically, we es- 
tablish a simulated sample of 28 sightlines from random 
sampling of the main sample, allowing duplications of 
individual sightlines. Then, we evaluate the cumulative 
iV(HI) distribution of the simulated sample. We repeat 
the procedure 10,000 times to determine the 68% scatter 
of T around the mean value in each 7V(HI) bin. 

3. THE METHOD 

The optical depth of Lyman limit photons along in- 
dividual lines of sight is determined according to tll = 
ctll xiV(HI), where <tll = 6.28 x 10~ 18 cm 2 is the photo- 
ionization cross section of hydrogen atoms. In principle, 
considering a sample of random sightlines from the star- 
forming regions in a galaxy together yields an estimate 
of the optical depth averaged over all viewing angles. In 
practice, we consider an ensemble of random sightlines 
toward GRBs in distant star- forming galaxies. The mean 
escape fraction of Lyman limit photons averaged over all 
directions is evaluated according to 

i—n 

(/e S c) = -Vexp[-(TLL^(HI)], (1) 

n ' 
i=l 

where the sum extends over the total number of n GRB 
sightlines in the sample. For our main sample presented 
in Figure 1, we find n = 28 and Equation (1) yields 
(/esc) = 0.02 ± 0.02. The error is estimated using the 
bootstrap re-sampling method described in § 2 and rep- 
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resents the 68% uncertainty in the mean value. We also 
determine a 95% c.l. upper limit (/ eS c) < 0.075. 

Parameterizing the cumulative A (HI) distribution by 
T[< JV(HI)] = A[N{}il)/N Q ] a , we find based on a X 2 
analysis log A = -0.58 ± 0.05 and a = 0.32 ± 0.03 for 
log iV = 20.5 over log iV(HI) = 16.5-21.5 (dash-dotted 
line in Figure 1). Equation (1) is expressed as 

r°° dT 
(/esc)=/ dNHi—-j- — exp[— ctll iV(H I)] 

JO flJVHI 

We determine (/ esc ) = 0.020 ± 0.003, where the errors 
represent the 68% uncertainties. 

4. DISCUSSION 

We have applied the ./V(HI) measured from the Lya 
(and in some cases Lyman series as well) absorption 
strength in early-time afterglow spectra to constrain the 
mean escape fraction in distant star-forming galaxies. 
Different from conventional methods to search for trans- 
mitted Lyman continuum photons, our estimated (/ eS c) 
is not subject to background subtraction uncertainties 
and does not depend on the intrinsic ultraviolet spectral 
shape of the host galaxies or dust distribution in the host 
ISM. In addition, it does not depend on the stochastic 
IGM Lya absorption along the sightlines toward these 
GRBs. Here we discuss possible biases in our estimated 
(/esc) due to the selection of GRB sightlines and impli- 
cations of our result. 

4.1. Observational Biases 

Our analysis considers only galaxies that host a GRB 
event. The presence of a GRB indicates that the ISM 
immediately surrounding the burst is photo-ionized by 
the bright afterglow. The intensity of the afterglow ion- 
izing radiation is capable of photo-ionizing all H I gas to 
a distance of r = 10 to 30 pc (Drain & Hao 2002; Watson 
et al. 2007; Prochaska et al. 2007b) that can exceed the 
typical Stromgrcn radius of an H II region. It is conceiv- 
able that the presence of a GRB may reduce tll along 
the sightline and the constraint on (/ esc ) can be consid- 
ered as an upper limit. 

In addition, our constraint is derived explicitly at the 
912-A Lyman limit transition. atoms. Additional absorp- 
tion due to dust would further reduce the estimated 
(/esc)) but is almost negligible at frequencies beyond 1 
Ryd (see Gnedin et al. 2007). However, some fraction of 
the GRB afterglows are optically "dark" and missed in 
our sample because afterglow spectroscopy was impossi- 
ble. While the majority of these are presumably associ- 
ated with high-z events (where the IGM absorbs most of 
the optical photons) or intrinsically faint afterglows, the 
remainder are associated with highly dust-extinguished 
sightlines (e.g. Pellizza et al. 2006; Rol et al. 2007). 
These extremely dusty sightlines are optically thick to 
ionizing radiation further strengthening the (f csc )- 

Intrinsically faint UV afterglows could occur in a par- 
ticularly low-density environment (Kumar & Panaitescu 
2000) and, in turn, preferentially probe optically thin 
sightlines. We note, however, that there is no notable 
correlation between the observed iV(HI) and the after- 
glow UV luminosity. For example, GRB 021004 and 



GRB 050820 have very different AT(HI) but comparable 
afterglow UV luminosity (e.g. Lazzati et al. 2006; Ves- 
trand et al. 2006). Furthermore, the isotropic equivalent 
energies of the GRB sample (L; so in Table 1) display 
a wide dispersion at all redshift, supporting that these 
GRBs do not preferentially originate in low-density en- 
vironment (e.g. Piran 1999). We argue that the pre- 
selection of bright UV afterglows has not restricted the 
analysis to a special sub-sample of GRBs. 

4.2. Implications and Future Work 

The association between GRBs and massive stars that 
typically have a short lifetime indicates that the GRB 
events occur close to the locations where their progeni- 
tor stars were formed and directly trace the current star 
formation rate. Late-time imaging surveys to search for 
the host galaxies of long-duration GRBs have shown that 
these GRBs occur in sub-L*, galaxies (Le Floc'h et al. 
2003; Sollerman et al. 2005; Fruchter et al. 2006) that 
exhibit on average large specific star-formation rates, 
SFR/Lb, (Christensen et al. 2004). This can be un- 
derstood as low-mass galaxies undergoing early genera- 
tions of star formation (Erb et al. 2006). The derivation 
of / osc from GRB sightlines reveals the escape fraction 
during the lifetime of the most intense ionizing sources, 
and applies to 'normal', star- forming galaxies that dom- 
inates the cosmic UV luminosity density. The low-mass 
nature also allows a direct comparison with predictions 
from high-resolution cosmological simulations (the dot- 
ted curve in Figure 1; Gnedin et al. 2007). 

Our 95% c.l. upper limit of (/ eS c) is comparable to 
the low values previously reported based on observa- 
tions of Lyman continuum photons from more luminous 
star-forming galaxies at z ~ 3 (Giallongo et al. 2002; 
Fernandez-Soto et al. 2003; Shapley et al. 2006). Addi- 
tional absorption due to dust would further reduce the 
estimated (f CS c), but is almost negligible at frequencies 
beyond 1 Ryd (see Gnedin et al. 2007). We estimate the 
total contribution of ionizing photons from sub-L* galax- 
ies, adopting the UV luminosity function determined for 
luminous star-forming galaxies at z ~ 3 from Adelberger 
& Steidel (2000). We derive a co- moving luminosity 
density at 1500 A of 2.2 x 10 26 h erg s" 1 Hz" 1 Mpc~ 3 
for sub-L* (0.1 — 1 L*) galaxies 7 . Applying an extinc- 
tion correction to the observed 1500-A flux (the authors 
estimated / CSC (1500A) « 0.2) and assuming an intrin- 
sic flux ratio between rest-frame 1500 A and 900 A of 
/(1500)//(900) = 3 adopted by Steidel et al. 2001 (but 
see Siana et al. 2007), we estimate a co-moving emissivity 
at 1 Ryd from sub-L* galaxies of < 2.8 x 10 25 h erg s^ 1 
Hz" 1 Mpc~ 3 for the 95% c.l. upper limit (/ esc ) < 0.075. 
To estimate the QSO contribution to the ultraviolet 
background radiation, we adopt the z = 3 QSO lumi- 
nosity function estimated by Hopkins et al. (2007). We 
find that the contribution to the ionizing background 
from QSOs of bolometric luminosity Lboi > 10 8 L Q is 
« 5 x 10 24 h erg s _1 Hz -1 Mpc~ 3 . While the uncer- 
tainties in these various numbers are large, this exercise 
shows that QSOs and sub-L* galaxies with (/ eS c) = 1 — 2 
% can contribute a comparable amount of ionizing pho- 
tons to the ultraviolet background radiation at z ~ 3. 

7 We adopt a A cosmology, S!m = 0.3 and S7a = 0.7, with a 
dimcnsionless Hubble constant h = _H"o/(100 km s _1 Mpc -1 ). 
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A larger sample of GRB sightlines with known A(H I) 
measurements for the host ISM is needed for improving 
the uncertainties in (/ esc )- In addition, follow-up imaging 
surveys for unvailing the emission properties of the GRB 
host galaxies are valuable for testing model predictions of 
star formation at high rcdshift. Specifically, cosmological 
simulations show that low-mass galaxies are inefficient 
in emitting ionizing radiation (Gncdin ct al. 2007). The 
luminosity distribution of these GRB hosts offers a direct 
test of these models. 

Finally, it is worth exploring whether a similar analysis 
could be performed with core-collapse supernovae (SN). 
These events also trace the death of massive stars over a 
broader mass range and provide a bright probe of the op- 
tical depth through the host galaxy. While the extreme 
line-blanketing in the far-ultraviolet of SN spectra may 
preclude the study of Hydrogen absorption, it may be 
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plausible to pursue an analysis of ISM metal-absorption, 
e.g. via the Mgll doublet. Indeed, the HST spectrum of 
the Type II SN 1999cm shows narrow absorption lines of 
Fell, Mgll, and Mgl (Baron et al. 2000). These features 
imply that NGC 1637 is optically thick to ionizing radi- 
ation along this particular sightline. A survey of z > 0.2 
SN could be carried out with a modest resolution, blue- 
sensitive spectrometer on a 10 m-class telescope. 
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TABLE 1 

The Sample of GRBs at z > 2 
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Note. — References: 1. Castro et al. (2000); 2. Castro ct al. (2003); 3. Vreeswijk et al. (2006); 4. Prochaska et al. in preparation; 5. 
Vreeswijk et al. (2004); 6. Jakobsson ct al. (2004); 7. Watson ct al. (2006);8. Berger et al. (2006b);9: Chen ct al. (2005); 10. Prochaska et 
al. (2007a); 11. Kawai et al. (2006); 12. Piranomonte et al. (2007); 13. Quimby et al. (2005); 14. Piranomonte et al. (2006); 15. Fynbo ct 
al. (2006b); 16. Cucchiara, Fox & Berger (2006); 17: Berger ct al. (2006a); 18. Price et al. (2007); 19. Cenko et al. (2006); 20. Berger & 
Gladdcrs (2006); 21. Savaglio et al. (2007); 22. Jakobsson et al. (2006); 23. Rol et al. (2006); 24. D'Elia et al. (2006); 25. Ruiz-Velasco et 
al. (2007); 26. Fynbo et al. (2006a); 27. Berger (2006); 28. Jaunsen et al. (2007); 29. Jakobsson et al. (2007); 30. Thoene et al. (2007a); 31. 
Berger et al. (2007); 32. Thocnc ct al. (2007b); 33. Hjorth et al. (2003); 34. Shin et al. (2007). 

a Isotropic equivalent energy release of 7-ray photons (erg s _1 ). For pre-Swift bursts, the energy interval corresponds to 200 — 2000 kcV 
(Bloom et al. 2003; Sakamoto ct al. 2005). For Swift bursts (2005 and on), the energy interval corresponds to 10 — 150 keV, and no 
k-corrections have been applied. 

b f a : flag for available Lya coverage in the afterglow spectra: = No coverage; 1 = Ly« coverage. 

c li: cumulative flag describing the properties of the ISM surrouding the GRB: 1 = presence of low ions (e.g. Sill, CII); 2 = the gas is 
optically thick at the Lyman limit; 4 = The Lyo absorption strength indicates JV(HI) > 2 X 10 20 cm~ 2 . For example, a flagigM = 7 
indicates the system is a DLA with associated low ions. 



